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ABSTRACT 

Internal gravity waves (IGWs) are naturally produced by convection in stellar 
envelopes, and they could be an important mechanism for transporting angular 
momentum in the radiative interiors of stars. Prior work has established that 
they could operate over a short enough time scale to explain the internal solar 
rotation as a function of depth. We demonstrate that the natural action of IGWs 
is to produce large scale oscillations in the solar rotation as a function of depth, 
which is in marked contrast to the nearly uniform rotation in the outer radiative 
envelope of the Sun. An additional angular momentum transport mechanism is 
therefore required, and neither molecular nor shear-induced turbulent viscosity 
is sufficient to smooth out the profile. Magnetic processes, such as the Tayler- 
Spruit dynamo, could flatten the rotation profile. We therefore conclude that 
IGWs must operate in conjunction with magnetic angular momentum transport 
processes if they operate at all. Furthermore, both classes of mechanisms must be 
inhibited to some degree by mean molecular weight gradients in order to explain 
the recent evidence for a rapidly rotating embedded core in the Sun. 

Subject headings: stars: interiors — Sun: rotation — waves 



1. Introduction 

During their pre-main sequence contraction, young solar-type stars are spun up to rota- 
tional velocities of the order of 100 kms" 1 . However, during their subsequent main sequence 



department of Astronomy, The Ohio State University, 4055 McPherson Laboratory, 140 West 18th 
Avenue, Columbus, OH 43210; dpa@astronomy.ohio-state.edu, pinsono@astronomy.ohio-state.edu. 

2 0n leave from Sobolev Astronomical Institute of St. Petersburg State University, Universitetsky Pr. 28, 
Petrodvorets, 198504 St. Petersburg, Russia. 

3 High Altitude Observatory, National Center for Atmospheric Research, P.O. Box 3000, Boulder, CO 
80307-3000; kmac@hao.ucar.edu. 



-2- 



(MS) evolution the surfac e rotation slows down as a result of an gular momentum loss through 



magnetized stellar winds (lKawalerlll988l ; lMatt fc Pudritzll2008l ). If a convective envelope and 



a radiative core of a solar-type MS star rotated independently of one another then the surface 
spindown would lead to a strong differential rotation be neath the envelope. In contradiction 



with this, helioseismic data (e.g.. ICouvidat et al.ll2003l ) reveal that the solar radiative core 



to the radius r 



of young cluster stars ( 


Stauffer & Hartmann 


Krishnamurthi et al. 


1997 


:lBarnesll2003 


) re- 



quires that the internal differential rotation only persists for timescales of the order of 20 
Myr to 100 Myr. These data indicate that, in radiative interiors of solar-type MS stars, 
there is an efficient mechanism of angular momentum redistribution that couples their core 
and envelope rotation. Unfortunately, its physical nature remains elusive in spite of many 
attempts to understand it made in the last years. Matters are further complicated by the 
requirement that an appropriate model of angular momentum transport in the solar-type MS 
stars should also reproduce the intricate variations (depletion) of the su rface Li abundance 



as fu nctions of age and effective temperature observed in the same stars (ISestito fc Randich 



2005|). 



A brea kthrough in solving this com plex problem has recently been announced by lTalon fc Charbonnel 



(120051 ) and I Charbonnel fc Talon! (120051 ) . They have shown how internal gravity waves (here- 
after, IGWs) generated by the envelope convection can extract angular momentum from 
rapidly rotating radiative cores of the solar-type MS stars on short enough timescales to 
explain both the spindown of young cluster stars and the quasi-solid-body rotation of the 
Sun. Furthermore, they have found that the quick flattening of the internal rotation profile 
by IGWs reduces element diffusion coefficients associated with hydrodynamic instabilities 
induced by differential rotation to values consistent with those constrained by the Li data. 

Pursuing the g oal of uncovering intrinsic causes of canonical extra mixing in low-mass red 
giant branch stars (jDenissenkov &: VandenBergi 120031) and trying to under stand the origin 
of fast rotation of red horizontal branch stars (jSills fc Pinsonneaultl |2000| ). we have been 
undertaking a critical review of different transport mechanisms in stellar radiative zones. 
When applying our test model of angular momentum redistribution by IGWs to a model of 
the present-day Sun, we have found that IGWs strongly disturb the solar internal rotation 
making it disagree with the helioseismic data. The disturbance can be eliminated only if 
there is another transport mechanism com peting with IGWs. I t s effic iency should exceed 



that of the rotational shear mixing used by iTalon fc Charbonnell (120051 ) by more than three 
orders of magnitude. Alternatively, a magnetic transport mechanism could be a competitor 
for IGWs. But in that case the mechanism preventing IGWs from disturbing the solid-body 
rotation of the Sun could itself be responsible for both shaping the internal rotation of the 



- 3- 



solar-type MS stars and assisting in the depletion of their surface Li abundance. This paper 
presents a discussion of computational results that support our conclusions. 



2. Angular Momentum Transport by IGWs 



In a gravitational field, any perturbation exerted on a fluid excites in it both acoustic (p- 
modes) and internal gravity waves (g- modes). The driving force for the latter is the buoyancy, 
as opposed to the pressure for the p-modes . A ge neral discussion of their combined linear 
theory can be found in the book by lLighthilll (119781 ). Our application of the IGW theory first 
addresses the following three questions: How do IGWs propagate through radiative layers? 
How do they redistribute angular momentum in a rotating star? How are IGWs generated? 
Answering the last question necessarily includes a discussion of the IGW spectrum. We start 
general comments on IGWs, then attempt to examine the underlying issues. 

The physics of angular momentum transport and chem i cal m i xing by IGWs in stellar 



radiat iv e zones has b een c o mprehen s ively de scribed bylPresa (11981 uGarcfa Lopez fc Spruit 




ZahnetaL (1997) 



Kim fc MacGregod (|200lh 



Rineot 



(19 m. Kumar et all (1199911 



Talon et al 



Montalban fc Schatzman 



(I2002h . and iTalon fc Charbonnell (120051 ). In 



all of these papers, IGWs are considered to be generated by large-scale turbulent fluid mo- 
tions either in the convective envelope or at the interface between the radiative core and 
convective envelope. The net energy flux of IGWs at the core/envelope interface r c is usu- 
ally estimated as Fe{t c ) ~ A4 t F c , where F c = p c v^ ~ 0.1 (L/4nr^) is the convective flux, 
L is the star's luminosity, and Ai t = u c /N c is the turbulent Mach number. In the last 
ratio, N c is the buoyancy frequency immediately beneath the interface, while u c = 2ttv c / A c 
is the turnover frequency of the largest convective eddy approaching the interface with the 
velocity v c . Everywhere in this paper u denotes the circular frequency. Although it should 
be measured in rads^ 1 , we will always express its values in units of /zHz, assuming that 
1 /iHz= 10~ 6 rads -1 . In the mixing- length theory (MLT) of convection, that we apply, the 
diameter of the largest turbulent eddy, which also measures its mean free path A c , is an «mlt 
fra ction of the local pressure scale height Hp. We employ the stellar evolution code described 
by lDenissenkov et al.l (120061 ). Gravitational settling is not included because it does not affect 
the IGW propagation. Indeed, IGWs can easily penetrate inner radiative cores of low-mass 
MS stars (ITalon fc Charbonnell 120051 ) where radial variations of the mean molecular weight 
fi are m uch stronger than those incurred from the operation of gravitational settling. We 
use the Grevesse fc Noels Jl993h mixture of heavy elements. Our calibrated solar model 
reproduces the solar luminosity (L & = 3.85 x 10 33 ergs _1 ) and radius (R & = 6.96 x 10 10 cm) 
at the solar age of 4.57 Gyr; this procedure yields the helium and heavy-element mass frac- 
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tions Y = 0.273, Z = 0.018, and the mixing length ckmlt of 1.75. Our solar model has 
r c « 0.713^0 and A^ t ,® ~ 8.3 x 10~ 4 , so that the energy luminosity of IGWs at the bottom 
of its convective envelope L E (r c ) = 47rr 2 F E (r c ) comprises only ~ 0.0083% of L . 

On their way from the bottom of the convective envelope toward the center of a solar- 
type MS star, IGWs experience radiative damping. This can be taken into account by 
applying a wave attenuation factor exp(— r) to L E (r c ). To calculate the effective optical 
depth, we use the relation 



T 



[1(1 + 1)] 3/2 



K 



dr' 



(a/N) 2 r' 3 



derived bv IZahn et all (11997I ). Here, N 2 

9* r 



N^ + N 2 , where 



N 2 



ad 



rad } 



and Nf. 



dr 



are the T- and //-component of the square of the buoyancy frequency. In the last two expres- 
sions, V r ad and V a d are the radiative and adiabatic temperature gradients (logarithmic and 
with respect to pressure), and g is the local gravity. The quantities 6 = — (<91np/<91nT) p 
and <f = (<91np/<91n/i) pr are determined by the equation of state. In our IGW computa- 
tions, the equation of state for the ideal gas is used. In this particular case, 5 = <p = 1. Also 
in equation (CQ), K = 4acT 3 /2>np 2 Cp is the radiative diffusivity with k and Cp representing 
the Rosseland mean opacity and the specific heat at constant pressure, respectively. The 
quantities / and a are introduced below. 

Turbulent eddies with different length and overturn time scales present in the convective 
envelope can generate a whole spectrum SE(r c ,l,m,u) of IGWs with different spherical 
degrees I > 1, azimuthal numbers m (\m\ < I), and frequencies u. If the star rotates and its 
angular velocity Q varies with r then the optical depth ([T]) depends on all three of the wave's 
spectral characteristics, the latter two entering it through the doppler-shifted frequency 
a = oj — m[Q(r) — 0(r c )] (we use a cartesian coordinate system in which the z-axis is colinear 
with the vector ft). Values of this frequency should be watched to remain between and 
N. When a — > 0, the optical depth r approaches the infinity, and the wave is completely 
absorbed by th e surrounding medium. On the other hand, when a — > N, the wave is totally 
reflected back (lRingotlll998l ). If only the IGWs with 1 < I < Z max and cu m i n < u < u; max are 
excited and propagate into the radiative core then their net energy flux at the core/envelope 
interface is 



F E (r c 



EE 

1=1 m=—l 



S E {r c , I, m, u)duj = M t F c . 



(2) 
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Besides energy, IGWs can also carry angular momentum the flux of which is 

Fj = -F E (3) 
a 



(IRingotl Il998l ; iKumar et al.lll999l ). Since both F E and a are positive, this relation means 
that prograde waves (those with m > 0) carry positive angular momentum while retrograde 
waves (m < 0) transport negative momentum. Note that before Ringot's elucidating paper 
the improper (negative) sign was used in equation ([3]) by many researchers creating some 
confusion in the field. 

Let us choose a frame of reference co-rotating with the convective envelope and denote 
Afi(r) = tt(r) — 0(r c ). We will assume that, in spite of the action of the Coriolis force in this 
noninertial frame of reference, the spectrum of IGWs generated by the envelope convection 
is still axisymmetric, i.e. S E (r c ,l,m,u) = SE(r c ,l, —m,u); if necessary, corrections du e to 



the Coriolis force can be taken into account later on (e.g., iTalon fe Charbonnell 120031 ). If 
AQ(r) = 0, i.e. if the whole star rotates uniformly, then a = lj, and r does not depend on 
m. In this case, the total angular momentum luminosity associated with IGWs 

Lj{r) = 4irr 2 V" V" m / S E (r c , I, m, u) exp{-r(r, I, m, u)} — (4) 
' * ' * /, , to 

1=1 m=-l J "n>in 

is equal to zero at any radius r below the convective envelope. 

On the other hand, if AO(r) 7^ 0, i.e. if the star rotates differentially, then IGWs will 
experience selective damping in the radiative core. Indeed, let us assume for example that 
Afl(r) increases with depth (when r decreases). In this case, the frequency a of a prograde 
wave with spectral characteristics I, m, and u will be doppler shifted to values smaller than 
the intrinsic frequency u, while that of a retrograde wave with the characteristics /, — m 
and u will become larger than u as the waves propagate inward. Hence, the difference 
between their attenuation factors exp{— r(r, I, m, u)} — exp{— r(r, I, — m, u)} will diminish 
with the depth. This happens because the prograde waves carrying the positive angular 
momentum experience stronger damping (if AQ(r) > 0) on their way into the radiative core 
than the retrograde waves carrying the negative angular momentum. The positive angular 
momentum absorbed by the surroundings at the beginning of the waves' path will spin up 
rotation locally compared to Q(r c ) + AQ while the negative momentum accumulating in 
the waves as they advance inward will be deposited deeper where the retrograde waves get 
eventually absorbed. This briefly sketched physics of angular momentum redistribution by 
IGWs is incorporated into the following PDE: 

2 dAtt 1 d ( 4 dAtt\ 3 1 dLj 



^ dt r 2 dr dr J 87r r 2 dr ' ^ 
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where v is a viscosity. Supplemented with eqs. ([I]), (j3J), an expression for the IGW spectrum 
at r = r c , and appropriate initial and boundary conditions, this equation describes how the 
star's internal rotation profile evolves in the presence of IGWs and a viscous force. 



3. Spectra of IGWs 



Given that the optical depth r, that determines the efficiency of damping of an IGW, 
strongly depends on the wave's spectral characteristics I, m, and uj, knowing the spectral 
energy distribution for IGWs at the core/envelope interface is as important as estimating 
their net energy. In this paper, we will use two analytical prescriptions (eqs. M and [7] below) 
for the spectra of IGWs. They correspond to two different physical mechanisms of the IGW 
excitation by turbulence in the convective envelope. The formal lower and upper limits for 
the intrinsic frequency of these IGWs are ui mi-n = oj c and u> max = N c . 

Unfortunately, multidimensional hydrodynamic simulations of turbulent convection in 
the solar envelope and its penetration into the radiative core give contradictory results on 
the spectru m of IGWs generated i n the se numerical experiments. For example, in their 2D 
simulations iRogers fc Glatzmaierl (120051 ) have found an IGW energy flux evenly distributed 
in frequency at least for I < 20 with a peak energy that is t hree orders of magni tude smaller 
than that predicted by equation (j7j). On the other hand, iKiraga et al.l (120051 ) claim that 
the broad frequency IGW spectrum reported in their earlier publication was an artefact 
of the 2D approximation. They emphasize that the IGW flux obtained in their new 3D 
simulations is of the same ord er as the one calculated from the simple parametric model of 
Garcia Lopez fc Spruitl (1199 lh based on the MLT. This finding encourages us to use spectrum 
dSJ) as our primary IGW model. 



3.1. The Garcia Lopez & Spruit Spectrum 



This prescription approximates the spectrum of IGWs generated at the core/envelope 
interface as a result of dynamic hitting at the radiative side of the interf ace by break - 
ing convective eddies. This appro ximation was prop osed and developed by iPressi ( 1l98ll ). 
Garcia Lopez fc Spruit <jl99lh . and IZahn et al.l (119971 ) . We use an expression derived in the 
last work 



S E (r c ,l,m,uj) = -p c v'l 



1 

2? 



(6) 



2 r ^ c N c l c oo L 

Equation ([6]) is obtained under the following assumptions: (i) the dynamic pressure in the 
waves of a frequency uj matches that produced by the convective eddies with the overturn 
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time ~u; , (wj the kinetic energy spectrum of the convective motions in the envelope is 
represented by the Kolmogorov law; (Hi) besides waves at their own length scale A, the con- 
vective eddies also excite IGWs with horizontal wave lengths Ah > A by the superposition of 
their incoherent acti on on the interface; these long er waves have a velocity amplitude reduced 
by the factor A/ Ah ( iGarcfa Lopez fc Spruitl Il99lh . In this prescription, / max = l c (to / tu c ) 3 ^ 2 , 
where l c = 2iir c /\ c . Applying the integration and summation ([2]) to the spectrum ([6]) and 
noticing that u c <C N c , we find that F E (r c ) = M. t F c as expected. 



3.2. The Goldreich et al. Spectrum 



The second prescription originated from the investigation of the stochastic excit ation 
of p-modes by turbulent convection in the Sun carried out by iGoldreich et al.l (11994J ) . In 
this model, IGWs are generated by the fluctuating Reynolds stresses pr oduced by turbulent 
fluid motions in the convective envelope. Following iTalon et al.l (120021 ). the emerging IGW 
spectrum at the core/envelope interface can be estimated as 



S E {r c ,l,m,uj) 



1^! f Rc a £_ 

2/47T.L r 2 
x exp 



dr J ^ ^ \ dr 



3\4 



V 3 \ 



l + (^r A ) 15 / 2 ' 



(7) 



where R c is the outer radius of convective envelope, T\ = X/v is the overturn time of con- 
vective elements of the size A = auurHp moving with the velocity v (at r = r c , A and v 
coincide with the parameters A c and v c defined in §[2]), and h u = Amin{l, (2lut\)~ 3 / 2 }. Like 
equation ([6]), the latter equation has been derived under the assumption that the turbulent 
motions in the convective envelope obey the Kolmogorov law. 

The radial displacement wave function £ r and the horizontal one 1 (1 + l)£h ; the latter 



being related to the former by the continuity equation (IZahn et al.lll997l ). are normalized to 



( Press 


1981; 


Zahn et al. 


1997; 


Kumar et al. 


1999) 



e(r)=C[/(/ + l)] 1 / 4 - 



exp 



-f 



\N\ y/ 1(1+1) 



dr' 



uj(pr\N\y/ 2 



(8) 



This equation shows that in the convective envelope, where iV 2 < 0, IGWs are evanescent. 
Therefore only those of them that are excited close to the core/envelope interface will effec- 
tively contribute to the IGW flux at r = r c . The constant C is adjusted in our computations 
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for the spectrum © to yield F E (r c ) = M t F c with F c = 0.1 (L/Airrl). The factor (2/)" 1 
in eqs. (JSJ) and ([Tj) takes into account the assumed energy equipartition between the wave 
counterparts with opposite signs of the azimuthal number (— I < m < I). 

In Fig. [TJ we have plotted the logarithms of the IGW energy luminosity summed over 
all available values of m at the core/envelope interface in our solar model. Those were 
calculated using the spectra (jHJ) (solid curve) and flTj) (dashed curves) with the same value 
of A4 t = 8.3 x 1CT 4 that gives L E (r c ) fa 3.2 x 10 29 ergs -1 . Apparently, they have quite 
different dependences on the spherical degree and frequency. Whereas the first spectrum 
does not depend on I at al£3 (eq. Ej), the second is estimated to be proportional to l' p , where 
p fa 1.6 for the first 4 degrees. Besides, the first spectrum declines with increasing uj much 
slower (with a power —3) than the second spectrum (with a power ~ —4.5). As a result, in 
the second case the dipole wave (1 = 1), that experiences the least damping (eq. [Tj), carries 
much less energy toward the radiative core than it does in the first case, especially at higher 
frequencies (compare the lower dashed and solid curve). 



3.3. Uncertainties in the MLT 

Fig. [2] shows how the convective and buoyancy frequency, u c and N, vary with the radius 
on the opposite sides of the core/envelope interface in our solar model. From this figure, it is 
evident that neither the minimum frequency of IGWs u m i n = u c nor the ratio Ait = u c /N c 
estimating the net energy flux Fe(t c ) fa M.tF & (dotted curve in the figure shows that the 
latter approximation becomes valid at r > r c + 0.5Hp), can be predicted with a confidence 
by the MLT of convection employed by us. The uncertainties are caused by the rapid growth 
of lo c and iV with an increasing distance from the interface. Strictly speaking, both uj c oc v c 
and N oc | V ra d — V a d| 1//2 should vanish at r = r c . However, as soon as we step aside from the 
interface, both of them jump up to finite values. In stellar model computations, their first 
nonzero values depend on spatial resolution: the higher the resolution is, the smaller these 
values are. For example, in our computations of the present-day Sun's model, the results of 
which are plotted in Fig. [2j we find oo c fa 0.30 /iHz and iV c fa 360 /iHz, hence M t ~ 8.3 x 10~ 4 . 
If we had taken into account convective overshooting beyond the formal lower boundary of 
convective envelope located at the radius r = r c , where V ra d = V a d (the Schwartzschild 
criterion), then we would have obtained a larger value of N c . In a value of A4 t = u c /N c , 
the increase of iV c due to the overshooting can partly be compensated by choosing a larger 
representative value of uj c that should account of its rapid growth with radius on a length 



The factor (21) 1 disappears after the summation over all m = — /, . . . , I. 
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scale much less than Hp immediately above the core/envelope interface (dashed curve) 



Besides the aforementioned uncertainties in the choice of representative values for uj c and 
N c , the MLT does not account for the fact, established by laborato ry experiments, observa- 



tions in the Earth's atmosphere and numerical simulations (e.g., seelStein fe Nord lund 



1989 



Cattaneo et al.fl~991 ; Rieutord fc Zahn 1995 ; Montalban fc Schatzmanl 200ol ; Rogers fc Glatzmaier 



20061 ). that downward flows in a strongly stratified convection zone are much more energetic 
and confined in their horizontal extent than upward flows. This may result in an underesti- 
mate of_thel_G^Vn^t_energy flux, especially in stars with deep convective envelope, like the 
Sun fairaga et al.ll2005h . 



Given the uncertainties in the MLT and the fact that we want to investigate the stability 
of the solar uniform rotation against perturbations by the IGWs that were supposed to be 
powerful enough to produce that uniform rotation in the past, we tentatively choose u min = 
u c = 0.30 /iHz while considering the turbulent Mach number A4 t as a free parameter with 
its minimum value equal to our estimated ratio uj c /N c rs 8.3 x 10~ 4 . This may actually put 
a conservative lower l imit o n the energetics of such IGWs because other authors, including 
Charbonnel fc Talonl (120051 ) who demonstrated how the IGWs could shape the Sun's uniform 
rotation, employed larger values of Ai t and sometimes also a higher u> m i n . Therefore, their 
low-degree high-frequency waves, that can penetrate deep into the radiative core, carried 
more kinetic energy than they do in our basic case, assuming, of course, that we use the 
same IGW spectrum. For the convective parameter l c , we use an MLT value l c = 30 from 
our calibrated solar model. 



4. Qualitative Description of Expected Solutions 



Aft er the proper (positive instead of negative) sign in relation (j3J) had been defined by 



Ringotl (119981 ) it became clear that the addition of IGWs to other processes responsible for the 
redistribution of angular momentum in stellar radiative zones would not be a trivial problem. 
Indeed, instead of resulting in an exponen tial decay of perturbations of the internal rotation 
profile toward a solid body rotation (e.g., iKumar fc Quataertlll997l ). the damping of IGWs 
leads to progressively growing deviations of local rotation from its initial state, even if this 
state is close to the solid body rotation. This is caused by the fact that in a region rotating 
faster than the convective envelope the prograde waves experience stronger damping than 
the retrograde waves, which results in a deposit of positive angular momentum there and 
hence in a further spin-up of this region. The opposite is true for a region rotating slower 
than the convective envelope. Of course, the perturbations of the rotation profile cannot 
grow infinitely large. First of all, an excess of (positive or negative) angular momentum in 
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them can slowly be dissipated through the molecular viscosity. Second and more important, 
when the rotational shear produced by the perturbations becomes strong enough to trigger a 
shear instability the turbulent viscous friction associated with the shear-induced turbulence 
will start to contribute to smoothing the perturbations out. 

Apparently, the outcome of the competition between the disturbing action of IGWs 
and the smoothing effect of viscous friction depends on their relative strength. Let us take 
a constant viscosity v = u . It is obvious that, as long as u is kept extremely large, the 
rotation profile will remain stationary because any perturbation of Q by IGWs will be quickly 
neutralized by the viscous dissipation. Of course, this does not preclude gradual temporal 
changes of the rotation profile as a whole due to the redistribution of angular momentum by 
IGWs as described by eq. fl5j). 

If we begin to decrease z/ then at some critical value of it, which is proportional to 
the local IGW energy flux, regular oscillations of Q will set in. In the Earth's atmosphere, 
this bifurcation of a stationary solution toward an oscillatory solution is believed to be 
observed experimentally as the quasi-biennial oscillations (QBO) of the mean zonal wind 
in the equatorial stratosphere disturbed by IGWs coming d own from the troposphere (e.g., 
Lindzen &: Holtonlll968l ; |Plumblll977l ; lYoden fc Hortonlll988l ) . In models of the solar-type MS 

bund as a solution 



stars, a behavior of the f2-profile resembling that of the QBO h a s been 
of equ a tion ([5]) nea r the t o p of radiative core by Kumar et al.1 (|19 99); 



(2001) 



Talon et al.l J2002h; Tralon fc Charbonnell (I2003h . and Icharbonnel fe Talonl (j2005h . 



Kim &: MacGregor 



Talon fc Charbonnell (120051 ) have called it the shear layer oscillations (SLOs). They have 
proposed that the SLOs work as a filter for IGWs on their way toward the radiative core. 
If Q increases with the depth in the core, as is expected in a solar-type MS star losing 
its angular momentum from the surface through a magnetized stellar wind, then the SLOs 
predominantly filter out the retrograde waves. These are absorbed closer to the center. 
Possessing the minimum angular momentum, it is the very central part of the core that 
slows down the first. Follow ing it, increasingly more an d more distant from the center layers 
get successively decelerated ( ITalon fc Charbonnelll2005l ). Note, however, that this theoretical 
prediction has recently been challenged b y results of a p r elimin ary analysis of the GOLF data 
on solar g-mode oscillations reported by I Garcia et al.l (120071 ). They suggest that the solar 
inner core at r < 0.15 R & may rotate three to five times as fast as the rest of the radiative 
zone. If these results are confirmed by a further analysis they will likely rule out the angular 
momentum redistribution in the Sun by IGWs. 

If we continue to reduce z/ further on then the oscillat ions will grow up in amp litude 



until finally they will turn into chaotic variations of Q (e.g., iKim fc MacGregor! 1200 ll ). 
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5. Results: Two Spectra, Five Viscosities 

In this section, we present and discuss results of our numerical solutions of the PDE (jSJ) 
that have been obtained with the IGW spectra (jSJ) and ([7]) using the viscosity prescriptions 
summarized in Appendix B. Besides restricting the azimuthal number m by the values of ±Z, 
we have also considered a limited set of the spherical degree. In most of our computations, 
we have only used a set of numbers I — 1, 2, 3, and 4. We do not think that adding higher 
degrees would qualitatively change our results and conclusions. Indeed, the optical depth in 
the wave attenuation factor exp(— r) increases as r oc I 3 for I > 1, therefore at a same radius 
in the Sun's radiative core a contribution to the local energy flux of a wave with a higher I is 
made from a part of the spectrum at a higher frequency to ~ Z 3 / 4 (solid curves in Fig. [3]). As 
both of the IGW spectra quickly decline with a growth of lu, our neglect of waves with the 
higher spherical degrees is unlikely to lead to a serious mistake. The same argument justifies 
our choice of the limited frequency interval 0.3 /zHz < uj < 3 /iHz because waves with higher 
frequencies transport negligible amounts of kinetic energy and angular momentum. 

Unlike the other papers cited in the preceding section, in this work we have not tried 
to solve the full problem of shaping the uniform rotati on of the solar radia t ive co re by 



IGWs. A tentative solution of it has been provided by ICharbonnel fc Talon! (120051 ). We 
look at this problem from another perspective. Let us assume that the redistribution of 
angular momentum in the Sun's core has alread y established its solid body outer envelope 



rotation as revealed by the helioseismic data (e.g.. lCouvidat et al.ll2003l ). Despite this, IGWs 
still continue to be generated by the envelope convection. We cannot expect a priori that 
the stability of the fi-profile against perturbations by IGWs strongly depends on its shape. 
Therefore, we consider it worth investigating whether these perturbations are sufficiently 
weak in the present-day Sun for them not to disturb noticeably the Sun's uniform rotation, 
provided that the IGWs producing these perturbations had enough energy in the past to 
couple the Sun's core and envelope rotation. Given our reformulation of the problem, we 
have chosen the following initial and boundary conditions for eq. ([5]): 

AQ(0,r) = Tc - T Af](0,r b ) for r b < r < r c , (9) 
r c ~ r h 

and AVI (t, r c ) = 0, Afi (t, r b ) = Afi (0, r b ), (10) 

where r b = r c — 0.6 R & . Thus, we assume that AQ initially increases with the depth linearly 
up to its maximum value Af2(0,r b ). Taking into account that for the present-day Sun 
Q (r c ) « 2.9 /iHz, the initial rotational shear in our computations decreases with the depth 
z = (r c — r)/R G as q = (<91nf2/<91nr) rs 3.1(x c — z)u max , where x c = r c /R Q , and M max = 
AQ (0,r b )/10~ 6 . In our basic parameter set, we will use a value of M max = 10~ 4 /iHz^C fi(r c ). 

To solve eq. ([5]) one needs to know stellar structure parameters, such as p, N, K, and 
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others, as functions of radius and time. However, because the internal structure of the Sun 
has not changed appreciably in the last billion years, we will use our model of the present- 
day Sun as a background for all of our IGW computations and we will watch that the total 
integration time in each of them does not exceed ~ 1 Gyr. We have solved eq. (jSJ) using an 
original method described in Appendix A. 



5.1. Constant Viscosity 

Fig. H] shows our results obtained for different values of the constant viscosity z/o,n = 
u/10 n using the spectra ([HD (panel a) and ([7]) (panel b). This figure illustrates the afore- 
mentioned bifurcation from stationary to oscillatory solutions that occurs at u 08 Ri5x 1CT 4 
in panel a and at z/q,8 ~ 8 x 10~ 4 in panel b, a sequence of oscillatory solutions (these are 



analogs of the SLOs discussed by lTalon fc Charbonnelll2005l ) for 5 x 10 5 < z/q,8 < 5 x 10 /: 



in panel a and for 7 x 10~ 5 < uq S < 8 x 10~ 4 in panel b, and a transition to the chaotic 
behavior of AQ at u 8 = 3.8 x 10~ 5 in panel a and at = 4.9 x 10~ 5 in panel b. Although 
in these computations we have only taken into account IGWs with frequencies from the 
narrower interval 0.3 /xHz < uj < 0.6 /iHz, this truncation does not depreciate our results 
because waves with these low frequencies carry 75% and 91% of the total energy for the 
spectral distributions ([6]) and ([7]), respectively. A comparison of panels a and b shows that 
the use of either of the two IGW spectra leads to similar qualitative results. 

Since the low-frequency waves get absorbed very close to the core/envelope interface, 
we have taken the zooming parameter k = 2.6 (see Appendix A) in order to resolve the short 
lengthscale variations of the rotation profile produced by them. For this k and for the used 
value of n = 8, we get a time scaling factor of 9.7yr for our dimensionless PDE (1A3|) . The 
real timescales of the oscillations of AQ are somewhat longer than this because, additionally, 
they are inversely proportional to the amplitude of the second term on the right-hand side 
of the PDE which is of order 10 -2 -10™ 3 . Therefore, the minimum time intervals between 
consecutive curves in Fig. H] are 10 3 yr (panel a) and 10 2 yr (panel b). A decrease of vq 
leads to both longer timescales and larger amplitudes of the oscillations because, in order to 
compete with the disturbing action of IGWs, a smaller viscosity needs a stronger shear to 
be built up, which means longer viscous dissipation times. 

So far, we have used our basic parameter set that includes a limited number of / = 1, 2, 3, 
and 4, and uses w max = 10 _4 /iHz in the initial and boundary conditions fl9l- [T0|) . Panel a in 
Fig. [5] demonstrates that neither the addition of 4 extra I values nor the increase of u max by 
the factor 10 3 , which results in the initial shear q m 0.22 near the interface, change much 
the period and amplitude of the oscillations of AQ. 
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Although the changes caused by the increase of Afi(0, n,) (w m ax) turn out to be unim- 
portant for our investigation of the ability of IGWs to disturb the i ntern al uniform rotation 
of the present-day Sun, they were shown by iTalon &: Charbonnell (120051 ) to be a matter of 
great importance in the problem of angular momentum extraction by IGWs from the radia- 
tive core of a young solar-type MS star. Indeed, a steep initial Af2-profile results in the SLOs 
that are asymmetric with respect to the line AQ = (red curves in Fig. [5^,). Such SLOs 
may work as a filter that predominantly absorbs prograde waves. This means that among 
low-degree high-frequency waves that arrive at a rapidly rotating central part of the star 
retrograde waves transporting negative angular momentum will be over-represented. Hence, 
when being damped in the core they will spin it down. 



5.2. Shear-Induced Viscosity 

The constant viscosity has been adjusted by hand to get one of the three possible out- 
comes of the competition between the disturbing action of IGWs and the smoothing effect 
of the viscous force on the rotation profile. It turns out that the quasi-periodic oscilla- 
tions of AQ are settled only when vq takes on a value from a rather narrow interval. It 
is unlikely that this accidentally happens in real stars. For the IGW filter composed of 
quasi-regular temporal and radial variati ons of AQ near the top of radiative core to work as 



proposed by ITalon fc Charbonnell (120051 ) some self- regulating mechanism for adjusting the 
proper viscosity values should apparently be operating there. One such mechanism could 
be a shear-induced viscous friction. In this case, a shear is readily built up as a result of 
selective damping of IGWs. Absorbed waves deposit their angular momentum locally, thus 
pushing rotation away from its stationary state. The viscosity coefficient (1B2[) used by us 
i s appropriate for descri bing mixing due to the rotation-induced secular shear instability 



(IMaeder fc Meynetlll996l ). It develops more easily than the dynamical shear instability but 
it acts on a smaller length scale It, such that a turbulent eddy of size ~/ t can effectively 
exchange heat with its surroundings while it travels a mean free path of order ~ l t . 



Talon fc Charbonnell (120051 ) have implemented this mechanism as follows. They aver- 
aged the turbulent diffusion coefficient (1B6I) over a complete oscillation cycle as well as over 
a radial extent of the SLOs using a Gaussian of width 0.2Hp. Thus obtained stationary 
viscosity profile was then used in their IGW computations. As is said in Appendix B, we 
actually use the same viscosity (lB2p but we allow it to vary with time and we do not average 
it over radius. For the same choice of IGW parameters as in the preceding section but for 
the shear-induced viscosity (IB2I) calculated with f v = 1, results of our solution of the PDE 
(IA3I) with the spectrum ([6]) are plotted in Fig. [5b. 
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It is importa nt to note tha t, like a few other pub l icatio ns (e.g.. lKim fc MacGregorll2001 



Talon et al 



2002h . the work of iTalon fc Charbonnell (120051 ) only contains a discussion of the 



SLOs near the core/envelope interface. Unlike them, we have decided to address the question 
whether the SLOs die out at greater depth or not. In order to shorten our computation time 
(i.e. in order to allow longer time steps) when solving the PDE (1A3|) in the bulk of radiative 
core we have bounded the IGW frequency by the values 0.6 /^Hz < uj < 3 /iHz, i.e. we 
have cut off the IGWs that produce the SLOs very close to the interface, like those shown 
in Fig. \5jp. In spite of th is, our minimum frequency still approximately equals the lowest 
frequency 0.5 /iHz used by lCharbonnel fc Talon! (120051 ). It is important to note that the IGW 
spectra we use are still normalized by equation (T5]) over the whole frequency and spherical 
degree intervals: 0.3 /iHz < uj < 3 ^Hz, and 1 < I < l c — 30. Because of its rapid decline 
with an increase of uj a spectrum normalized with a; m i n > 0.3 /iHz would have more energetic 
high-frequency waves than ours, hence it would produce even stronger oscillations of Afl 
than those obtained by us. 



We have solved equation (1A3I) in the k = 1-zoomed depth interval < z < z^ = 10 k 



X 



(x c — xt,), where x c = 0.713 and x^ = x c — 0.6, for time periods less than 1 Gyr. Hoping 
to relate, later on, the angular momentum transport and element mixing by IGWs in the 
Sun to canonical extra mixing in low-mass red giants, we have taken the parameter / v = 
20 in the expression (1B2I) for the shear-induced viscosity because with about that value 



Denissenkov et al.l (120061 ) succeeded in reproducing evolutionary abundance variations of Li 



and C in the atmospheres of cluster and field red giants. We have found that, even after 
having been enhanced by this large factor, the shear-induced viscosity fails to extinguish 
large scale SLOs deeper in the solar radiative core. The blue curve in Fig. [6^ represents 
an envelope of oscillation amplitudes of AQ. The farther inward from the core/envelope 
interface, the higher the AQ oscillation amplitude and the longer its characteristic time are. 
In the outer half of the radiative core the maximum amplitudes by far exceed the deviations 
of the helioseismic data (red squares with error bars) from the uniform rotation profile. A 
better agreement with the experimental data is obtained if we choose f v = 10 3 (purple curve 
in panel a). We have tried this value as well because the 3D hydrodynamic simulations 
by iBriiggen fc Hillebrandtl (120011 ) of turbulent mixing induced by the shear instability have 
shown that equation (IB1I) may underestimate the coefficient of turb ulent diffusion by three 



orders of magnitude. Of course, it is not clear whether the results of IBriiggen fc Hillebrandt 



(120011 ) can be applied directly to real stars, given that they were obtained assuming plane- 
parallel geometry and without taking into account "the effects of rotation, nuclear reactions, 
and variations in radiative processes" . 



However, before making any conclusions from the results of these computations we have 
to check if the viscosities induced by the shear flows outlined in panel a are high enough for 
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the turbulent fluid motions producing them not to be broken down by the molecular viscosity. 
For this to be true, the flow Reynolds num ber Re = v w J u mo i must exceed the critical Reynolds 
number Re c ~ 40 (jSchatzman et al.ll2000l ). In panel b, we have plotted a time averaged u v for 
the cases of f v = 20 (blue curve) and f v = 10 3 (purple curve). In the same plot, green curve 
depicts z/ mo i while red curve presents the quantity Re c x v mo i. Comparing blue, purple, and 
red curves in panel b, we conclude that the shear-induced turbulence can only be sustained 
near the base of the convection zone where both the shear q is sufficiently strong thanks to 
the very short lengthscales of the SLOs and the quantity proportional to z/ v /^ 1 (fig) -2 (eq. 
BT]and dashed curve in Fig. [3]) steeply increases with r. This raises the following important 
question: what alternative shear dissipating mechanism works in the Sun's outer radiative 
core that successfully competes (as follows from the helioseismic data) with the disturbing 
action of IGWs? 



5.3. Molecular Viscosity and Ohmic Diffusivity 



Fig. [6] shows that the viscosity needed to counteract the distortion of the solar rotation 
profile by IGWs should not necessarily be too high. Taking into account that the blue and 
purple curves in Fig. [6b represent the time averaged u v , whose real values change with time 
following the oscillations of Afi, it seems worth testing if the molecular viscosity can smooth 
out the large scale SLOs alone. The shear-induced turbulent viscosity can only be used 
down to a depth z ~ 0.04-0.05 because below this region the ratio ^v/^moi becomes smaller 
than the critical Reynolds number (Fig. [6b). It is interesting that the size of this region 
approximately coincides with the thickness of the solar tachocline in which the latitudinal 
differential rotation of the conne ctive envelope is transformed into the quasi-soli d body ro- 
tation of the radiative core (e.g., ISpiegel fc Zahnlll992l ; ICharbonneau et al.lll998l ). We want 
to find out if the viscous force in our computations fails to reduce amplitudes of the SLOs 
in the outer radiative core to values consistent with the helioseismic data simply because we 
use the IGW spectrum ([6]) instead of ([7]). 

It is possible that this negative result is a function of the assumed IGW spectrum. In 
order to respond to these questions, we have employed the spectrum ([7]) and a combined 
viscosity v = z/ v (/v)+/moiX^moi, where z/ v (/ v ) is a substitute of equation (1B2I) in which we set 
/ v = 20 for < z < 0.04 and f v = for 0.04 < z < zy,- We have als o investigated models with 
more vigorous IGWs. Furthermore, as lCharbonnel fc Talon! (120051 ) claimed, the total energy 
luminosity of IGWs produced by fluctuating Reynolds stresses in the convective envelope of 
their solar model was 8.5 x 10 29 ergs -1 . This is about 2.7 times as large as our estimated 
value of Le{t c ) ~ 3.2 x 10 29 ergs -1 . Therefore, we have increased our turbulent Mach number 
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A4 t = 8.3 x 10~ 4 by this factor and renormalized the spectrum ([7]) respectively. Results of 
these computations obtained with the factor / mo i = 2 are plotted with purple curves in 
Fig. [7J For test purposes, we have also repeated these computations with an extended set of 
the spherical degree I — 1, 2, . . . , 7, 8 (blue curves). Fig. [7] shows that the molecular viscosity 
(even after it has been dou bled) cannot compete with the IGWs that have been shown by 



Charbonnel fc Talon! (120051 ) to be powerful enough to shape the Sun's solid body rotation. 



On the other hand, it turns out that a viscosity proportional to the ohmic diffusivity 
v — /mag x ?7mag can extinguish the IGW-induced SLOs everywhere in the solar radia- 
tive core except the tachocline region for a value of / mag > 10 (Fig. [8]). This seemingly 
pure academic exercise has some sense. Let us assume that differential rotation in the so- 
lar radiative core has been s uppre s sed by magnetic process es, e.g. lik e thos e proposed by 



Charbonneau h MacGregorl (119931 ) , iMenou fc Le Merl (120061 ) , or ISpruitl (119991 ) (for more de- 



pi 
r( 



tails on the latter, see next section). However, in order that magnetic fields generated by 
these processes not to decay too quickly through the ohmic dissipation, an effective magnetic 
diffusivity r\ c associated with them must exceed ?7 ma g. 



5.4. Effective Magnetic Viscosity 



Spruitl (119991 ) has proposed a magnetohydrodynamic mode of angular momentum trans- 
port in radiative zones of differentially rotating stars. Fluid elements experience large-scale 
horizontal displacements caused by an unstable configuration of the toroidal magnetic field 
(one consisting of stacks of loops concentric with the rotation axis). Small-scale vertical 
displacements of fluid elements are coupled to the horizontal motions, which can cause both 
mild mixing and much more effective angular momentum transport. Spruit's key idea is 
that no initial toroidal magnetic field is actually needed to drive the instability and mixing 
because the unstable field configuration can be generated and maintained by differential 
rotation in a process similar to convective dynamo. The Spruit dynamo cycle consists of 
two consecutive steps: first, a poloidal field is generated by the vertical displacements of 
the unstable toroidal field; second, the new poloidal field is stretched into a toroidal field by 
differential rotation. 



The Spruit mechanism produces a huge magnetic viscosity u e > 10 9 cm 2 s 1 ( iDenissenkov fc Pinsonneai 



20071 ) that could indeed prevent IGWs from disturbing the solar uniform rot ation. However 



it cou ld produce that uniform rotation itself, without being assisted by IGWs (jEggenberger et al 



20051 ). 



The original prescription for the effective magnetic diffusivity and viscosity in the model 
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of S pruit's dynamo h as recently been criticized by iDenissenkov fc Pinsonneaultl (120071 ) (see 
also Zahn et al. 2007 ). The principal critical argument is that Spruitl ( 1999 ) has overesti- 
mated the horizontal length scale of the Tayler instability that causes the concentric mag- 
netic loops to slip sideways. Spruit assumed the length was of order a local stellar radius. 
Denissenkov fc Pinsonneaultl (120071 ) suggested that one has to account of the Coriolis force 
when estimating the instability's horizontal length scale. This reduces the diffusivity by 
about three orders of magnitude and results in the following expression: 



K o 2 

^ ~ 2 Tat \2 1 ' cm s • 
l iv rj_3 



(11) 



There is no need to do any further computations to understand that, with this revised dif- 
fusivity, magnetic fields generated by the Spruit dynamo in the solar radiative core would 
immediately be dissipated through the ohmic diffusivity, hence the whole transport mecha- 
nism would not function. Indeed, equation (fTTj) gives r/ e « 1.8 x 10 2 cm 2 s -1 < r/ mag near the 
base of the solar convection zone even if we take q ~ 1 which obviously exceeds the upper 
limit constrained by the helioseismic data. The revised prescription may only work to reduce 
differential rotation in a model of the young Sun in which both f2_6 and q have much larger 
values (Denissenkov et al., in preparation). 



6. Conclusion 



Our numerical solutions of the angular momentum transport equation (|5j or IA3j) have 
demonstrated that neither the molecular viscosity nor the shear-induced turbulent viscosity 
can reduce the large scale oscillations of angular velocity in the solar outer radiative core 
caused by selective damping of IGWs, provided that the net energy flux of these waves is 
strong enough to shape the Sun's solid body rotation. If waves were the sole mechanism, we 
would therefore expect to see large deviations from rigid rotation. The amplitudes of these 
oscillations are found to be too large to agree with the helioseismic data. Our result holds 
even when the molecular and shear-induced viscosity are multiplied by large factors. We have 
proved that only a viscosity exceeding the ohmic diffusivity by a factor of > 10 can smooth out 
the IGW-induced oscillations of the rotation profile. This may be an indirect indication that 
some magnetic processes are at work here. To be more precise, our finding actually satisfies 
a necessary condition for such processes to work because the effective magnetic diffusivity 
r] e associated with them must exceed the ohmic diffusivity 77 mag . Otherwise magnetic fields 
generated by them will decay through the ohmic resistivity too quickly. For example, we 
have found that magnetic torques are strong enough to successfully compete with the action 
of IGWs only if the original prescription for the Taylor-Spruit dynamo, in which r] e 3> r/ ma g, 



-18- 



is used. However, this particular case turns out to be irrelevant to our p roblem because the 



Tayle r- Spruit mechanism can shape the solar solid body rotation alone (jEggenberger et al. 
20051 ). without being assisted by IGWs. 



The helioseismic data suggest that either there is an efficient angular momentum trans- 
port mechanism in addition to IGWs that smooths out the SLOs produced by the waves or 
the spectral energy distribution of IGWs is different (lower) from those used by us. Although 
the latter assumption leads us outside the scope of our formulated problem we will comment 
on it. It is possible that strong toroidal magnetic fields in the solar tach ocline filter out the 



IGWs with the doppler-s hifted frequency a above the Alfven frequency (IKumar et al.lll999 



Kim &: MacGregorll2003l ). For the minimum frequency 0.6 /iHz used in our computations of 
the large scale SLOs in the solar outer radiative core the magnetic field strength requi red to 
prevent the inward wave propagation is about (3 x 10 5 )// Gauss (IKumar et al.lll999l ). For 
1 = 1, this corresponds to quite a strong field. If it is present in the tachocline then the 
waves with high spherical degrees will be trapped there. However, we have only considered 
IGWs with uj > 0.6 /iHz and low 1 values. Most of them are likely to propagate below the 
tachocline. In order to find out if an enhanced viscosity in the tachocline can hinder the 
propagation of low-degree high-frequency waves into the solar radiative core we have done 
test computations in which the shear- induced viscosity was increased by a factor of 10 3 . 
Their results plotted with black curves in Fig. [7] show that this does not help to solve the 
problem. 

Another, more radical possibility is that the form of the IGW spectra employed by us 
is completely wrong. For instance, the IGW spectrum estimated in the 2D hydrodynamic 
simulations by iRogers fc Glatzmaierl (120051 ) has a flat energy distribution which goes three 
orders of magnitude below the peak luminosity in our Fig. [TJ Apparently, if we applied that 
spectrum in our computations then even the molecular viscosity could easily smooth out the 
SLOs produced by such IGWs. However, it is evident as well that IGWs with this energy 
distribution could not produce the uniform ro tation of the Sun by its pre sent age (multiply 
the ages of the rotation profiles in Fig. 1 from ICharbonnel fc Talon! 120051 by a thousand). 



To summarize, we do not see what microscopic or pure hydrodynamic processes could 
smooth out the large scale SLOs induced by IGWs in the solar oute r radiative core. There- 
fore, we agree with the conclusion made by lGough fc Mcintyrd (119981 ) about "the inevitability 
of a magnetic field in the Sun's radiative interior" . Indeed, if IGWs are as strong as described 
by our employed spectra then this magnetic field is required to trigger magnetic processes 
that will counteract the disturbing action of IGWs on the solar rotation profile. On the 
contrary, if IGWs are weak then we are in need of such magnetic processes to extract an 
excess angular momentum from the solar interior on the early MS. 
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During the preparation of this work, results of a preliminary analys is of the GOLF data 
on solar g-mode oscillations have been published by I Garcia et al.1 (120071 ) suggesting that the 
solar core at r < 0.15 R & may rotate three to five times as fast as the rest of the radiative 
core. If these results prove to be correct, they will seem to rule out the angular momentum 
redistribution in the Sun by IGWs, as proposed by ICharbonnel fc Talon! (120051 ). because in 
that case it would have been the Sun's inner core to be spun down first. If these results are 
confirmed, it will mean that the strong /x-gradient in the Sun's central region has prevented 
any angular momentum transport mechanism from operating there. There would also be 
st rong implicat ions for magnetic angular momentum transport, ruling out prescriptions (such 



as 



Spruit 



20021 ) that predict a weak sensivity to \i gradients. 
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A. Solution of the Angular Momentum Transport Equation 

In order to shrink the problem's parameter space, we consider a contribution to the 
transport of angular momentum only from the waves possessing the maximum possible az- 
imuthal number \m\ = I, i.e. from those carrying the maximum angular momentum, both 
positive and negative, at given values of I and uj. T o solve the main PDE (jSJ), we use a 



numerical method based on the ideas implemented by ISaravananl (Il990l ) in his model of the 



wind flow in the Earth's stratosphere influenced by tropospheric IGWs. 

Taking into account the assumed axial symmetry of the IGW spectrum, i.e. that 
Se{t c , I, rn, uj) = Se(t c , I, —m, uj), and the fact that the optical depth in the wave at- 
tenuation factor only depends on even powers of the doppler-shifted frequency a (eq. [TJ, we 
recast the net angular momentum flux for m = ±Z as 

L (r) ' maX /"^max 

F ^ r ) = -f L T = 12 l / Sj(l,u,)exp{-T(r,l,l,u>)}dLJ, (Al) 



where 



Sj{I,u) = S e {tMH)u- X - (A2) 



In the interval — u c < uj < u> c , we set Sj(l,u) = 0. Note that Sj(l, —uj) = —Sj(l,u). 
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From the computational standpoint, we find it convenient to convert eq. (jSJ) to the 
following dimensionless form: 

du 1 d f 4 du\ _ a nik f L \ dF 

di ~ px± dz \ px Vn dz J px* \lJ dz ' [A6) 

where x = r/R Q , and 

i q3 ' max r Umai r& 

F(z, u) = / s A l > u ) ex P{" / G («'> 1 > U ~ lu ) d2 '} duJ - ( A4 ) 

In the last equation, G(z, I, uu — lu) = [1(1 + 1)] 3//2 7 _1 (z) (u — lu)~ 4 , where 7 _1 (5) = 
10 _fc x 0.2065 x KqN_ 3 (N t ) 2 _ 3 (x c — 10~ fc z)~ 3 with x c = r c /R Q . We neglect the square 
root a/1 — (u/N) 2 in the integrand's denominator in eq. because we will only consider 
w max N c in which case \a\ = \u — lu\ < 2u max N. 

When deriving eqs. (1A3HA4]) . we have normalized our basic variables, which are 
assumed to be initially expressed in cgs units, as follows: u = Afi/10~ 6 , v n = u/10 n , 
z = 10 k x (x c -x)J= (10 n+2k /R 2 Q ) t, K 6 = K/10 6 , and A^_ 3 = A^/10- 3 . The quantity k is 
a sort of zooming parameter. Taking k > allows us to look with the scrutiny at results of 
IGW damping taking place close to the base of convective envelope, where variations of u 
may occur on a very short lengthscale of order 10~ 3 - 10~ 2 R & . Having done these transfor- 
mations, it turns out that a n ^ = 1.362 x lO 8- ™^, z — corresponds to the core/envelope 
interface, and t is measured in units of 1.535 x io 14 ~ n_2fe yr. The factor 10 3 in eq. (IA4I) 
comes about from a combination of the factor 0.1 that estimates the ratio of the convec- 
tive flux to the total flux from the star and the factor 10 4 that represents the reciprocal to 
the normalization constant for the turbulent Mach number that has been included into the 
coefficient a n ^. Given that u is measured in pHz, a n ^ has additionally been multiplied by 
the factor 10 6 . This is necessary to do because the wave frequency that we actually use is 
uj = uu/10~ 6 , therefore, when substituted into eq. (1A4I) . relation (1A2|) should be taken in the 
form Sj = 10 6 x S E (w/10- 6 )" 1 . 

In equations (IA3HA41) . all integrals are replaced with series of trapezoids while all 
derivatives are approximated by finite differences. The zoomed depth interval < z < ib 
(here, ib = 10 fc (r c — r^)/ R Q ) is divided into M equal subintervals by M+l mesh points while 
the axisymmetric frequency intervals — uu m3iX < uu < —uu c and uu c < uu < uu max , representing 
the retrograde and prograde waves, respectively, are divided into L subintervals each. The 
resulting system of M + 1 nonlinear algebraic equations is linearized assuming that after 
every integration time step At the r atio max{| Auj/uA }jL <C 1. In order to speed up the 



computations, we follow the idea of iSaravananl (119901 ) to interpolate functions containing 
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G(£j, — luj) in £j, I, and the combination (tUi — luj) using initially prepared and stored 



tables 



B. Viscosity Prescriptions 



The viscosity v in eq. (jSJ) plays a very important role because, depending on its value, the 
viscous friction either succeeds or not in smoothing out oscillations of the rotation profile 
in the radiative core growing in response to the local deposit of angular momentum that 
accompanies the absorption of IGWs. In this work, we employ 5 different prescriptions for 
v as well as some of their combinations. These are a constant viscosity uq, the molecular 
viscosity v mo \ (it dominates over the radiative viscosity in the solar-type MS stars), a viscosity 
proportional to the magnetic (ohmic) diffusivity ?7 ma g, a viscosity u v associated with vertical 
turbulence produced by the secular shear instability induced by differential rotation, and 
an effective viscosit y v e related to magnetic torques generated by the Tayler-Spruit dynamo 
JSpruitlll999l . l2002h . 



For t he viscosity due to t he sh ear- induced vertical turbulence, we use the expression 
derived by iMaeder fc Meynetl (119961 ) multiplying it by a free parameter f v 



/A- 



8 K 
fv x -Ri c j^tt 2 q 2 , 



(Bl) 



where Ri c — j is the critical Richardson number, and the shear q = (d In fl/d In r). The 
parameter f v takes into acco unt the fact that, according to hydrodynamic simulations by 
Briiggen fc Hillebrandtl (120011 ). the original prescr iption may underestimate the viscosity by 
the factor of ~ 10 3 . Alternatively, ICanutol (120021 ) has supposed that Ri c should be at least 
four times as large as its classical value. After the same normalization used in Appendix A, 
we have 



V^v)n = fv x ~ 10 
5 



2k-n x 2 



(N 7 



(du\ 



(B2) 



To compute the effective magnetic viscosity, we use both original Sprui t's equations (see 
our § !5.4p and the equations revised by iDenissenkov fc Pinsonneaultl (120071 ) 



r 2 VLrjl 



1/3 



(B3) 



2 We take advantage of the fact that |o;, — luj\ < 2ui n 
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where the effective magnetic diffusivity is 



K 

l U q . 



(B4) 



Combining the last two equations and normalizing the variables, we find 



2.686 x 10 5 x 103 



k—n 



x 



2/3 



fdu\ 



2/3 



(B5) 



It should be noted that we have implicitly assumed that N 2 = in eqs. (IBlj) and (1B4j) . 
This approximation may be valid in the outer part of radiative core if we neglect the [i- 
gradients produced by the gravitational settling and radiative levitation of chemical elements. 
Under this assumption, our choice of u v is equivalent to that made bvlTalon fc Charbonnel 
(120051 ) . Indeed, although they have used a prescription proposed bv lTalon fc Zahnl (1l997T l 



v v = -Ric—ril q 

5 Nj, 



;i + D h /K) 



l + (NZ/N*)(l + K/D h ) 



(B6) 



that takes into account a reduction of the stable thermal stratification in the radiative core by 
strong horizontal turbulence described with a diffusion coefficient D h 3> z> v , putting N 2 = 
in eq. (1B6I) and noticing that in evolve d solar-type MS stars D h K (e.g., see Fig. 14 in 
the paper of iTalon fc Charbonnel! 120051 ) transforms z> v into our u v . 
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Fig. 1. — Logarithms of the energy luminosity of IGWs at the Sun's core/envelope interface 
summed over all available azimuthal numbers (m = —I, ...,/) as functions of I and u. Solid 
curve corresponds to spectrum (jBj), dashed curves — to spectrum for I — 1, 2, 3, and 4 
(from the lower to upper curve). 
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Fig. 2. — The convective overturn frequency u c (dashed curve), the buoyancy frequency N 
(solid curve), and the ratio of the convective flux pv% to the total flux F Q (dotted curve) as 
functions of the relative depth (expressed in units of the pressure scale height) on the both 
sides of the Sun's core/envelope interface. 
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Fig. 3. — The frequency uj t =\ for which the effective optical depth (eq. [Q) in the IGW 
attenuation factor exp(— r) equals to one at a given radius in the case of uniform rotation 
(from the lower to upper solid curve: / increases from 1 to 4, and 0J T= x oc [/(/ + 1)] 3 / 8 ). Dashed 
curve is the normalized stellar structure parameter in expression (IBID for the shear-induced 
viscosity. 
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Fig. 4. — The shear-layer oscillations (SLOs) of AQ = Q(r) — Q c at the top of solar radiative 
core computed (eq. IA3I) with the constant viscosity z/ 08 = z^o/lO 8 using the IGW spectra 
(JHJ) (panel a) and ([7]) (panel b) for the following sets of parameters: 0.3 /iHz < uj < 0.6 /iHz; 
I = 1, 2, 3, 4; m — ±1; and M max = 10~ 4 /^Hz in the boundary and initial conditions (P-fTOl). 
Panel a: black curves — z/o,8 = 5 x 10~ 4 , the time interval between consecutive curves is 
At = 10 3 yr; blue — z/ 0i8 = 10" 4 , At = 5 x 10 3 yr; green — z/ 0j8 = 5 x 10" 5 , At = 5 x 10 3 yr; 
red — z/ ,8 = 3.8 x 10" 5 , At = 10 4 yr. Panel b: black — u 0>8 = 8 x 10~ 4 , At = 10 2 yr; blue 
- z/ 0>8 = 10" 4 , At = 10 3 yr; green — z/ 0j8 = 7 x 10~ 5 , At = 10 3 yr; red — z/ , 8 = 4.9 x 10" 5 , 
At = 2 x 10 3 yr. 
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Fig. 5. — The SLOs at the top of solar radiative core computed with the constant (panel a) 
and shear-induced (panel b) viscosity using the IGW spectrum (jBJ) and the same sets of basic 
parameters as in the previous figure. Exceptions are the enhanced value of M max = 0.1 fiHz 
for green and red curves, and the extended set of I — 1, 2, . . . , 7, 8 for red curve. In panel 
a, all curves have z/ 0i s = 10~ 4 , and At = 5 x 10 3 yr. Curves in panel b have been computed 
using the parameter / v = 1 in eq. ( ]B2j) . and At = 5 x 10 3 yr. 
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Fig. 6. — Panel a: blue and purple curves are the envelopes of amplitudes of the SLOs 
deep in the solar radiative core computed with the shear-induced viscosity (eq. IB2I with 
/ v = 20 — blue curve, and with / v = 10 3 — purple curve) using the IGW spectrum (jf)]) 
and the same sets of basic parameters as in Fig. HI except for the cut-off frequency interval 
0.6 /xHz < uj < 3 uHz. Red squares with error bars in panel a are the helioseismic data from 
Couvidat et al.l (120031 ). Panel b: green curve — the molecular vi scosity v m r,u red - — th e 



quantity Re c x u mo i, where the critical Reynolds number Re c = 40 (ISchatzman et al.l 120001 ); 
blue and purple curves — the time averaged shear-induced viscosities produced by the SLOs 
whose envelopes are plotted in panel a with the same color. 
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Fig. 7. — Same as in the previous figure but for spectrum (JTj) multiplied by the factor 2.7 
(see text) and for the combined viscosity v — z/ v (/ v ) + 2v mo i. Exceptions are blue curves 
that were computed for I = 1, 2,..., 7, 8. In the expression ( 1B2I) for 2/ v , the following 
parameters have been used: / v = 20 (blue and purple curves) and f v = 10 3 (black curve) for 
< (r — r c )/R Q < 0.04, while / v = (all curves) for (r — r c )/R Q > 0.04. For comparison, 
red curves are computed using spectrum (jHj). 
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Fig. 8. — Same as in the previous figure but for spectrum ([6]) and for the combined viscosity 
v = v v (/ v ) + / mag x ?7 ma g. In both cases, the parameter / v = 20 has been used in the 
expression (1B2j) for < (r — r c )/R & < 0.04. Blue curve is computed with / mag = 2, purple 
— with / mag = 10. Black curve shows the magnetic diffusivity r] mag . 



